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I. INTRODUCTION 

The long-lived doubly heavy baryons open a new bat- 
tlefield in the study of heavy quark dynamics in various 
aspects such as the interplay of strong and weak inter- 
actions in decays as well as the confinement (see recent 
review on this subject in [Q) . Indeed, there are some fea- 
tures of these baryons in contrast to the mixed flavor me- 
son B c composed by two heavy quarks, for which one can 
explore the technique of nonrelativistic QCD (NRQCD) 
§ in the form of potential NRQCD (pNRQCD) § or 
velocity-counting NRQCD (vNRQCD) §, |, |, while the 
description of doubly heavy baryons involves the inter- 
action with the light or strange quarks constituting the 
baryon together with the doubly heavy diquark. There- 
fore, one should combine the above NRQCD approaches 
with the heavy quark effective theory (HQET) || con- 
structed in order to take into account the interactions of 
local heavy spinor field with the soft light quarks and glu- 
ons. Thus, the doubly heavy baryons give the possibility 
to test these methods in different conditions specific for 
the baryons QQq or QQs. 

An experimental opportunity for such the studies could 
be quite real at Runll of Tevatron. Indeed, the B c me- 
son was observed by the CDF Collaboration at FNAL || 
under a lower statistics of Runl. Its spectroscopic char- 
acteristics and decays were theoretically described and 
predicted in the framework of QCD sum rules 10 , po- 
tential models [l2| and operator product expansion 
fl3f . Those predictions served to isolate some reasonable 
constraints preferable for the detection of B c signal such 
as the intervals of mass, production rate, lifetime and 
branching ratios of various decay modes with quite a high 



efficiency. The strategy on the study of doubly heavy 
baryons follows the same way as for B c fl4| . First, one 
estimated the spectroscopic characteristics in the frame- 
work of potential models 16 and sum rules of QCD 
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or NRQCD |L7|, that resulted in definite predictions of 
mass spectra and wave functions, which determine the 
normalization of cross-sections in terms of soft binding 
factor Jl8| . Second, corresponding yields in hadronic col- 
lisions and e + e~ annihilation can be obtained in the per- 
turbative theory (fourth order in a s for the production of 
two pairs of heavy quarks in QCD), so that the rates are 
large enough in order to expect a positive result of ex- 
perimental search for the doubly heavy baryons at Runll 
Pq . Third, the potential models and operator product 
expansion were used to calculate both exclusive and in- 
clusive branching ratios 20), respectively, as well as 
the lifetimes |^0| . A new fie d of interest for the inves- 
tigations is radiative, electromagnetic or hadronic, tran- 
sitions between the quasistable states in the families of 
baryons with two heavy quarks. A first step in the study 
of this problem was recently done in [ p2| , wherein some 
preliminary results were obtained on the electromagnetic 
transitions between the levels of Sf, c . 

In the present paper we analyze the basic spectroscopic 
characteristics for the families of doubly heavy strange 
baryons Qqq' = {QQ' s) in the framework of potential 
model and compare the results with the calculations in 
NRQCD sum rules and lattice simulations for the ground 
states. 

A general approach of potential models to calculate 
the masses of baryons containing two heavy quarks was 
in detail discussed in |l5|, |l6) . There are two clear phys- 
ical arguments in this problem, which differ by kinds of 
interquark forces in the baryon. The first motivation on 
the form of description is the pair interactions, while the 
second is the string-like presentation based on the Wilson 
loop for three static sources, where one source interacts 



2 



with other two sources through a string connected to the 
diquark. For the doubly heavy baryons, these two ap- 
proaches result in quite certain predictions, which can 
be clearly distinguished in the splitting between the ex- 
citations as described by Gershtein et al. in 16 . The 



peculiarity of doubly heavy diquark is its small size, that 
can be used in order to construct an effective approxi- 
mation in the form of combined NRQCD and HQET by 
implementing the hierarchy of following scales: the size of 
QQ'-diquark antitriplct color subsystem, tqqi , which is 
about inverse both the characteristic momentum transfer 
between the heavy quarks \k\ and the relative momen- 
tum of heavy quark motion \p\ ~ |fc| ~ mQ,Q' • v with 
a small nonrelativistic velocity v <C 1, the confinement 
scale h-QCD for the dynamics of light or strange quark 
and the heavy quark masses mg^i , so that 

r QQ' ■ Aqcd < 1, Aqcd < m Q . 

Under such conditions, we can factorize the doubly heavy 
diquark as a local source of colored field in the interac- 
tions with the light quark, while the dynamics in the 
diquark can be reliably described in terms of nonrela- 
tivistic spinors, i.e. NRQCD, pNRQCD or vNRQCD. In 
the present paper we explore the quark-diquark factoriza- 
tion of baryon wave-functions in the framework of non- 
relativistic potential model and use the QCD-motivated 
potential combining the known perturbative calculations 
for the static potential at short distances |2^, Q and the 
linear confining term at large distances in the way sim- 
ilar to the Buchmuler-Tye method [p5| extended to the 
three-loop (3 function for the effective charge in the static 
potential [GJ. Such the calculations differ from the ap- 
proaches of pNRQCD and vNRQCD, since the pNRQCD 
deals with the static potential of heavy sources fixed at 
a distance r not including the static energy of ultra- 
soft fields such as the contributions from the quark-gluon 
string or sea, while the vNRQCD describes the nonrela- 
tivistic coulomb system in the limit of negligible contri- 
butions by ultrasoft fields valid at ttiq^qi ■ v 2 3> Aqcd, 
that is broken for the heavy-heavy systems composed of 
bottom and charmed quarks under study. The static po- 
tential K 2 jj] is consistent with the high-virtuality nor- 
malization of coupling constant in QCD as well as with 
the slope of Regge trajectories, that define two scale pa- 
rameters of model. The heavy quarkonia spectra Q and 
leptonic constants |2(| are described with quite a good 
accuracy in this approach, that results in the fixed val- 
ues of heavy quark masses implemented in the potential 
model with the nonrelativistic quarks. The connection 
of these quantities with the pole and current masses of 
heavy quarks was discussed in [[I). 

We adjust the nonrelativistic Schrodinger equation for 
the description of light and strange quark dynamics. So, 
we determine the constituent mass of light and strange 
quarks in the way of minimal binding energy of quarks 
in the static field of colored source in the framework of 
presentation with the constituent mass formed as a piece 
of ultrasoft energy in the linear confining term of the 



potential or the string (see Section II. B). Then we test 
the estimated mass values of light and strange quarks in 
the potential model by calculating the masses of bound 
states with a single heavy quark, i.e. the charmed and 
beauty mesons. We discuss uncertainties and apply the 
procedure for the doubly heavy baryons. 

If the heavy quarks composing the antitriplct-color di- 
quark, have the same flavor, Q = Q', we have to take 
into account the Pauli principle for the identical fermions. 
Then we find that for the symmetric, spatial parity P- 
even wave functions of diquark, ^(r) (the orbital mo- 
mentum equals Ld — 2n, where n — 0, 1, 2 ... ) the spin 
wave-function of diquark should be symmetric too, and 
the summed spin of quarks equals 5 = 1, while for the 
antisymmetric, -odd functions 'l'd(r) (i.e. Ld = 2n + 1), 
we have S = 0. 

Under factorization of diquark and strange quark dy- 
namics, we accept the following notations for the clas- 
sification of levels in the system of QQs baryons: the 
summed spin of two heavy quarks, their orbital mo- 
mentum and principal quantum number are denoted by 
the capital letters S 1 , L and rid, respectively, so that 
the ground state of diquark, for instance, is marked as 
n 2 d +1 L — l 3 S, while the motion of strange quark is la- 
belled by a pair of lowercase letters n s I. Such the no- 
tations are based on the presentation that these quan- 
tum numbers are approximately conserved if we neglect 
the multipole interactions in QCD p7[ with the emission 
of soft gluons, which can be absorbed by the associated 
strange quark in the baryon or involved into the scatter- 
ing with the emission of kaon, for example, g + s — > q + K 
with q = u, d. In this expansion there are the following 
suppression factors mentioned above: the diquark size 
with respect to the scale of confinement, the relative mo- 
mentum of heavy quarks with respect to the heavy quark 
mass, that leads to suppression of both the chromoelec- 
tric and chromomagnetic dipole transitions, respectively. 
Moreover, in the doubly heavy baryons with identical 
flavors of two heavy quarks some transitions between 
the levels or mixing operators have a double suppression, 
since for the corresponding operators one could need the 
properties providing the change of summed heavy quark 
spin AS = 1 together with the change of their orbital 
momentum AL = 2m + 1, m G Z. Say, we expect that 
the excited 2 1 P level of diquarks bb and cc would be 
quasistable under the transition to the ground level l 3 ^. 
This picture can be not valid for the be diquark, since 
both values of summed heavy quark spin are admissible, 
and the spin changing operators are not removed and can 
result in a significant mixing of diquark states labelled by 
the summed spin. So, we restrict our consideration of Qb c 
mass spectrum by the spin-dependent splitting of ground 
state, while for the doubly heavy baryons fl cc and O^t. we 
present complete spectra of families. 

Next point is the following evident condition for the ap- 
plicability of diquark dynamics factorization. So, we cal- 
culate the size of diquark subsystem in order to demon- 
strate the reliability of calculations. We will see that 
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the low-lying levels have quite small sizes, especially in 
the case of bb diquark, while the higher excitations are 
large with respect to the distance to the strange quark, 
so that the approximation of local colored source for such 
diquark levels results in less reliable estimates of baryon 
masses. 

The paper is organized as follows. In Section II we 
describe the static potential explored in the paper and 
calculate the relevant constituent masses of light and 
strange quarks. We test the accepted approximations 
in estimates of masses for mesons with a single heavy 
quark. Then we introduce the spin-dependent forces. In 
Section III we present numerical results. Our conclusions 
are given in Section IV. 



II. NONRELATIVISTIC POTENTIAL MODEL 

Under the factorization of dynamics inside the doubly 
heavy diquark and in the system of strange quark and 
diquark, we use the nonrelativistic Schrodinger equation 
in order to solve the corresponding two-body problems, 
that yields spin- independent wave- functions and levels. 
We use the three-loop improved static potential in the 
Schrodinger equations with the heavy quark masses ad- 
justed to the observed data on the bottomonium and 
charmonium. The procedure for the determination of 
constituent masses of light and strange quarks with the 
given potential is described and tested with the data on 
the heavy-light mesons. Then we introduce the spin- 
dependent forces treated as perturbations. 



A. The static potential 

In QCD the static potential is defined in a manifestly 
gauge invariant way by means of the vacuum expectation 
value of a Wilson loop p8[ , 

V(r) = - lim iln(W r ), 



Wr = tr"Pexp ig 



(1) 



Here, T is taken as a rectangular loop with time exten- 
sion T and spatial extension r. The gauge fields are 
path-ordered along the loop, while the color trace is nor- 
malized according to tr(..) = tr(..)/trl. This definition 
corresponds to the calculation of effective action for the 
case of two external sources fixed at a distance r during 
an infinitely long time period T, so that the time-ordering 
coincides with the path-ordering. Moreover, the contri- 
bution into the effective action by the path parts, where 
the charges have been separated to the finite distance 
during a finite time, can be neglected in comparison with 
the infinitely growing term of V(r) ■ T. 

Generally, in the momentum space one rewrites the 
above definition of QCD potential of static quarks in 



terms of relevant quantity a v representing a so-called V 
scheme of QCD coupling constant as follows: 



V(q 2 ) 



-a 



4wa v {q 2 ) 



(2) 



After the introduction of o = the (3 function is 
actually defined by 
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so that in the perturbative limit at high q 2 the coefficients 
were calculated up to three-loop order in the V scheme 
by the two- loop matching of static potential with af s 
|| HJ. Then /32k = and 02 = ~ aiffi* + 
(02 — a 2 )/3Q ls , where the coefficients ai corresponds to the 
short-distance expansion 
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Note that expansion (||) cannot be straightforwardly ex- 
tended to higher orders of perturbative QCD because of 
infrared problems, that result in nonanalytic terms in the 
three-loop perturbative potential as was first discussed by 
Appelquist, Dine and Muzinich j28[ , 

Buchmiiller and Tye proposed the procedure for the 
reconstruction of f3 function in the whole region of charge 
variation by the known limits of asymptotic freedom to a 
given order in a s and confinement regime. Generalizing 
their method, the /3px function found in the framework 
of asymptotic perturbative theory (PT) to three loops, 
is transformed to the (3 function of effective charge as 
follows: 
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(5) 



where the exponential factor in the second term con- 
tributes to the next-to-next-to- leading order at a — > 0. 
This function has the essential peculiarity at a — ► 0, 
so that the expansion is the asymptotic series in o. At 
a — > oo the (3 function tends to the confinement limit 
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(6) 



which gives the confinement asymptotics for the static 
potential at long distances r as usually represented by 
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the linear potential (see discussion in ref . pgfl ) 

V coal (r) = k-r. (7) 

The construction of @ is based on the idea to remove 
the pole from the coupling constant at a finite energy, 
but in contrast to the "analytic" approach developed in 
p0| and modified in |Hj , the "smoothing" of peculiarity 
occurs in the logarithmic derivative of charge related with 
the (3 function, but in the expression for the charge itself 1 . 
Indeed, in the one-loop perturbation theory we have got 
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and the pole can be cancelled in the logarithmic deriva- 
tive itself, so that 
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As we see in the perturbative limit a — ► the deviation 
in the (3 function is exponentially small, and the usual 
solution for the running coupling constant is valid. 

Equation (||) can be integrated out, so that implicit 
representation of effective charge can be inverted by the 
iteration procedure, so that well approximated solution 
has the form 



a(p 2 



1 



A,ln(l + 7^ 2 )£ 



(8) 



where ??(/i 2 ) is expressed through the coefficients of per- 
turbative /3 function and parameter I in (|^), which is 
related to the slope of Regge trajectories a' P and the in- 
tegration constant, the scale A Q|. Thus, the dimension- 
less parameter / determines theoretically unknown ratio 
of perturbative scale in the QCD coupling constant A to 
the confinement scale involved by the Regge slope. 

The slope of Regge trajectories, determining the linear 
part of potential, is supposed equal to ot! p — 1.04 GeV~ 2 , 
so that in (0) we put the parameter k = 2 ^ a , . We use 

also the measured value of QCD coupling constant 
and pose 

af{m 2 z ) = 0.123, 

as the basic input of the potential. The transformation 
into the configuration space was done numerically in [Q , 
so that the potential is presented in the form of file in 
the notebook format of MATHEMATICA system. 



In |3l[ the analytic approach is extended to the j3 function, too. 



The analysis of quark masses and mass spectra of 
heavy quarkonia results in the following values ascribed 
to the potential approach |J : 



= 1.468 GeV, 



= 4.873 GeV. 



(9) 



Thus, the spectroscopic characteristics of systems com- 
posed of nonrelativistic heavy quarks are determined in 
the approach with the static potential described above. 



B. Constituent masses of light and strange quarks 

A fast moving light or strange quark interacting with 
a static source of gluon field nonperturbatively emits vir- 
tual quarks and gluons forming a quark-gluon sea or a 
string. Such the cloud of confined virtual fields with the 
valence quark is not a local object, which has some inter- 
nal excitations, of course. However, these excitations cor- 
respond to instable hybrid states related with additional 
qq pair or glueball degrees of freedom in the hadron. We 
do not consider such exotics in the present paper. Then, 
the fast moving light quark surrounded by the virtual 
soft and ultrasoft fields with no internal excitations can 
be described as a whole object with an effective mass, 
which we call a constituent mass [i q . Since we postulate 
the nonrelativistic Schrodinger equation, we fix the dis- 

2 

persion law of constituent quark: E q — 2—, Further we 
present a consistent determination of fi q , when the static 
potential entering the Schrodinger equation is fixed. 

The constructive procedure is the following. Since the 
constituent mass is a part of energy in the string confin- 
ing the quarks, and this confining energy is represented 
in the static potential by the term linearly growing with 
the quark separation, we argue that this part of energy 
should be subtracted from the potential. So, we will solve 
the Schrodinger equation for the light quark with the con- 
stituent mass n q and a static source with the potential 
V[r) — fio, where /zo is the energy subtraction. We ex- 
pect that the parameters fi q and /io are very close to 
each other, of course. A difference between them should 
be suppressed, since it is determined by such systematic 
reasons as the nonlocality of constituent quark as well as 
its dispersion relation. At fixed [1q we can investigate the 
dependence of binding energy of constituent quark on the 
parameter [i q . So, we numerically solve the Schrodinger 
equation 



P 

2n q 



V(r) 



*(r) = [A(/i,)+/io-M,Mr), (10) 



for the ground state. Thus, the binding energy is given 
by A(yU g ), and it is shown in Fig. El 

We see that the binding energy A has an optimal value 
corresponding to its minimum. We ascribe the position 
of minimum \i q = /i* = 0.37 GeV as the valid constituent 
mass of light quark attributed to the Schrodinger equa- 
tion with the given static potential. The difference be- 
tween the normalization point /j,q and the constituent 
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A GeV 




FIG. 1: The binding energy of light quark in the potential 
field of static source obtained from the solution of ( |l0| ) versus 
the constituent mass [i q . The normalization fio = 0.37 GeV 
corresponds to the minimum of binding energy. 



mass fi* can be extracted from the experimental data 
on the spin-averaged ground-state masses of heavy-light 
mesons by comparison with the theoretical expectation 



M Q (lS)=m Q + A(/x*)| 



<e!) 

2^Q 



+ 5/i, (11) 



where we put Sfi = /i* — /io an d use the subleading ap- 
proximation including the kinetic energy of heavy quark 
with a finite mass toq, whereas the average (p 2 ) can 
be estimated numerically with the wave functions ex- 
tracted from the Schrodinger equation. Remember, that 
(p 2 ) = 2/i* T, where T is the kinetic energy, which is phe- 
nomenologically independent of the quark flavor, since 
the bound states are posed in the intermediate region 
with the change of coulomb regime to the linear confine- 
ment, so that the potential is close to the logarithmic 
form, in which the kinetic energy is flavor-independent 
©■ 

The experimental data determine the spin-average 
masses in accordance with 



M Q (1S) 



3M V + M P 



where My p are the masses of vector and pseudoscalar 
states. Then, we find that (|ll|) is consistent with the 
experimental data if 



Sfi = 35 MeV, 



(12) 



and we see that Sfi -C fi* as expected. 

We remark that the same result on Sfi can be repro- 
duced by the numerical solution of relevant Schrodinger 
equation for the two-body problem with the correspond- 
ing kinetic terms for the heavy and light quarks, so that 



the difference between the values of Sfi in these two ap- 
proaches occurs less than 5 MeV, that points to the reli- 
ability of this result. 

The situation with the determination of constituent 
mass of strange quark is very similar, but slightly com- 
plicated. Indeed, the strange quark has a valuable cur- 
rent mass depending on the normalization point. For 
definiteness we fix this mass at a low virtuality rele- 
vant to the dynamics of bound states, so that we put 
m s = 0.24 GeV, that is consistent with the extraction of 
strange quark mass from the QCD sum rules at the scale 
1 GeV |34|] . Then we solve eq. ( |To| ) with the substitutions 
fi q — > fi s = m s + fi and po — > /xg, where fi should be close 
to /ig, and add the following perturbation: 



SV a = k ■ ((r) - (r )) 



(r) 



(13) 



where k is the coefficient in the linear term of poten- 
tial, the average size (r) is calculated under the wave- 
functions dependent of the scale fi s , so that (r) = 

1/2 

[/ ^^(r)^(r)r 2 d 3 r\ , and this perturbation has quite 
a clear origin. Indeed, the length of string having the 
weight /z s equals ^ s /k, and we ascribe this weight to the 
constituent mass. However, in the case of quark pos- 
sessing a valuable current mass, the hadron can have the 
size (r) , which is less than the attributed length of string 
fi s /k determining the constituent mass. Then, we have 
to subtract a part of fake energy that is given by the ex- 
cess of length in order to get a consistency of such the 
presentation. The value of (r*o) is determined at a nor- 
malization point /zg given by the minimum of binding 
energy for such the constituent strange quark. 
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FIG. 2: The binding energy of strange quark in the potential 
field of static source obtained from the solution of Schrodinger 
equation versus the constituent mass (is- The normalization 
— 0.29 GeV corresponds to the minimum of binding en- 
ergy. 

The dependence of binding energy in the field of static 
potential for the strange quark versus the constituent 
mass is presented in Fig. 0. 



6 



Thus, in the potential approach we put the position of 
minimum for the binding energy as the constituent mass 
of strange quark, so that 



IK 



= 0.53 GeV. 



(14) 



Following the same procedure as for the light quark 
we can use the experimental data on the heavy-strange 
mesons in order to extract the value of normalization 
point /ip entering as the difference Sfi s = fj,* — /Xg into 
the equation of ground state mass, so that using eq.(pd|) 
with (p 2 ) = 2T n*niQ I (fx* + tuq), we find 

5fi s = 15 MeV, (15) 

while the solution of Schrodinger equation taking into 
account the finite masses of heavy quarks in D s and B s 
mesons results in 



10 MeV, 



(16) 



which we will use in our numerical estimates. 

The above values of <5/i and 6fi s give the estimates of 
systematic accuracy of potential model under considera- 
tion, so that the uncertainty of our calculations for the 
masses of bound states is about SM = 40 MeV. 



k ' 



fill 
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mass is less than 0.42 GeV, then we deal with the condi- 
tions in the heavy-light meson, while at fi s greater than 
0.42 GeV, the situation with the strange quark takes 
place, since at low virtualities in the bound state the 
running mass of strange quark expanded in perturbative 
series over a s has a significant contribution of renormalon 
p5[ , reflecting the infrared singularity in a s . Then, the 
running mass of strange quark is greater than 0.42 GeV, 
and the correction ( |l3| ) is justified. The renormalon con- 
tribution can be subtracted from the running mass, and 
this can be done by a redefinition of S/i s in order to in- 
clude the renormalon contribution in the value ^Jq, as we 
have performed above, so that /ip is rather small. Thus, 
the running mass of strange quark m s = 0.24 represents 
the so-called subtracted running mass, that is correlated 
with the small value of /Zq, since we have rearranged the 
soft contribution between these two quantities. 

Next, the above constituent masses of light and strange 
quarks are fixed in the procedure for the ground states 
of mesons with the static heavy quark, and we do not 
perform the same procedure for the excitations and use 
the fixed values in the predictions for the doubly heavy 
baryons, not only for the ground states, but also for ex- 
citations, too. 

So, we numerically solve the following Schrodinger 
equations 



P . PQQ>' 
2^1* 2m 



Pq 2 , PQ' 2 



V(r) 



1 



V{r) 



tf 8 (r) = e s * s (r), (17) 



*QQ'( r ) = tQQ'^QQ'{ r ), 



2uiq 2itiqi 2 
where the mass of diquark is determined by 

m = rriQ + m Q > + e QQ > , 
and the baryon mass is equal to 

M = m + m s + Sfi s + e s . 



FIG. 3: The average size of heavy-strange meson versus the 
constituent mass /i s (the solid curve). The corresponding 
length of string with the weight /i s is also show by the dashed 
line. 



Thus, we completely determine the method for the cal- 
culation of spin-independent levels in the systems with 
the heavy, light and strange quarks. 



The difference between the conditions of constituent 
mass formation in the heavy-light and heavy-strange 
mesons is demonstrated in Fig. |3| where we show the 
dependence of average size of meson on the constituent 
mass in comparison with the length of string with the 
weight of constituent mass. We see that if the constituent 



C. Spin-dependent corrections 

Following ||, |3§], we introduce a specified form of 
spin-dependent corrections causing the splitting of nL- 
levels, so that in the system of heavy diquark containing 
the identical quarks we have 
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(18) 



3 3m^ 



Sq 2 4^(r) - \ ° 



3 m2 4L 2 , 



[6(£ d • S 1 ,,) 2 + 3(X rf • S d ) - 2i^2] 



where the first string corresponds to the relativistic cor- 
rection to the effective scalar exchange, while the second 
string represents the terms due to the single-gluon vector 
exchange with an effective coupling constant a s depend- 
ing on flavors of quarks composing the system. In ( 18 ) we 
take into account the color factor corresponding to the 
antitriplet state of diquark, i.e. we add 1/2 in front of 
usual expression for the quark-antiquark colorless state, 
and substitute the static potential V(r) for the color sin- 
glet sources. The last term in ( |l8|) represents the tensor 
forces expressed in terms of orbital and summed spin of 
quarks as was shown in fl2]| . 

Taking into account the interaction with the strange 
constituent quark can be done in an analogous way. So, 
we have to explore the following evident kinematics for 
the motion of two heavy quarks posed in the same point 
with the distance r to the strange quark: 

Sq ■ [r x p Q ] = ~S q - L, S Q -[rx p s ] = Sq ■ L, 



S s -[rx p Q ] 



1 



S s ■ L, S s ■ [r x pj = S s ■ L. 



In this kinematics, the first term appears in the interac- 
tion with the effective scalar exchange, the second stands 
in the exchange by the effective gluon field, and similar 
expressions appear in the terms with the spin of strange 
quark. 

Then we adopt the current-current form of interaction 
between the strange and heavy quarks with the appropri- 
ate antitriplet-state factor of 1/2 in front of static poten- 
tial and sum up the terms related with the heavy-strange 
subsystems, so that we can explore a usual technique for 
the derivation of spin-dependent perturbations similar to 
the Breit potential in QED, and for the interaction of S- 
wave diquark with the strange quark we get 




dV(r) 8 



rdr 



3 r 3 



O.s 



L S 



-S d -S s 4ir5(r) 



3 mem* AL 



[6(L • Sf 



3 m Q fi* r 3 3 3m Q/ uJ 
3(£ • S) - 2L 2 S 2 - 6(L ■ S d f - 3(i • S d ) + 2L 2 S 2 d ] 



(19) 



1 



where S = S d + S s . In ( p~9| ) we see that this form of 
spin-dependent forces coincides with the expression that 
could be derived under assumption of local doubly heavy 
diquark with the spin S d interacting with the strange 
quark, so that the perturbation in the quark-antiquark 
system with the diquark mass rriQQ = 2rriQ is exactly 
reproduced. 

The value of effective parameter a s can be determined 

by 



Us pvH2(T)m rcA /k 2 QCD y 1 1 

where /3q = 11 — 2nf/3 and rif — 3, m rc d is the reduced 
mass of quarks composing the two-particle system, T is 
the kinetic energy in the quark system, so that numeri- 
cally we get Aqcd ~ 113 MeV from the comparison of 



theoretical expression 

AM(ns) = l^^\R nS (0)\ 2 , (21) 

with the experimental data on the system of cc 

AM (IS, cc) = 117 ± 2 MeV, (22) 

where R n s(r) is the radial wave function of quarkonium, 
and it is calculated in the potential model under study. 

In the above estimates we explore the fact that the 
average kinetic energy of quarks in the bound state prac- 
tically does not depend on the flavors of quarks, and it 
is given by the following values 

(T d ) «0.19GeV, (23) 

and 

(T s > w 0.38 GeV, (24) 
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for the antitriplet and singlet color states, correspond- 
ingly. 

For the identical quarks inside the diquark, the scheme 
of LiS-coupling well known for the corrections in the 
heavy quarkonium, is applicable. Otherwise, for the in- 
teraction with the strange quark we use the scheme of 
jj-coupling. Then, L ■ S s is diagonal at the given J a , 
(J 8 = L + S s , J = J s + J), where J denotes the to- 
tal spin of baryon, and J is the total spin of diquark, 
J = Sd + La- 
in order to estimate various terms and mixings of 
states, we use the transformations of bases 



\J;J S ) = £|J;£)(-1)< J 



+S S +L+J) 



and 



x r>s -3)12./, \)\ J L S J j } . (I'M 



\J;J S ) = J2\J;J d )(-l)( J + s °+ L +V 
Jd 



V(2J d + l)(2J s + l)\ g j J f 



(26) 



where S = S s + J, and J d = L + J. 

For the hyper-fine spin-spin splitting in the system of 
quark-diquark, we use the presentation with the local di- 
quark for both the interaction of S'-wave diquark with 
the strange quark and that of P-wave diquark with the 
s-wave strange quark. Then we introduce the pertur- 
bation analogous to the spin-spin term in (fl9j) with the 
substitution of Sd — > J- 

Thus, we have defined the procedure of calculations for 
the mass spectra of doubly heavy baryons presented in 
the next Section. 



III. NUMERICAL RESULTS 

In this Section we calculate the mass spectra with 
account for the spin-dependent splitting of levels. As 
we have clarified in the Introduction, the doubly heavy 
baryons with the identical heavy quarks allow quite a reli- 
able interpretation in terms of diquark quantum numbers 
(the summed spin and the orbital momentum). Dealing 
with the excitations of 6c-diquark, we show the results on 
the spin-dependent splitting of the ground lS-state, since 
the emission of soft gluon breaks the simple classification 
of levels for the higher excitations of such diquark. 

The quark-diquark model of bound states is the most 
reliable for the system with the more heavy quark. 
Therefore, the calculations for S&h and Qbb are the most 
accurate, while the corrections due to the finite size of 
diquark can be valuable in S cc and fl cc . 

In what follows we consequently refer to the diquark 
and strange quark quantum numbers: a principal quan- 
tum number n and an orbital momentum, ridLdn s l s . The 



results on the characteristics of diquarks: the masses, 
sizes and wave functions, calculated with the potential 
K 2 [Q are very close to the values estimated in the 
Buchmiiller-Tye potential. The estimates are presented 
m Tables | and y. The wave functions and binding en- 
ergies of strange quark slightly depend of the diquark 



A. Qt,b baryons 

Denote the shift of level by A( J ) marked by the total 
spin of baryon J. So, for lS2p we have got 



= -10.5 MeV. 



(27) 



The states with the total spin J = | (or |), can have 
different values of J s , and, hence, they have a nonzero 
mixing, when we perform the calculations in the pertur- 
bation theory built over the states with the definite to- 
tal momentum J s of the strange constituent quark. For 



J = I , the mixing matrix equals 



-18.8 -4.7 



-4.7 39.4 



MeV, 



(28) 



so that the mixing practically can be neglected, and the 
level shifts are determined by the values 



A' ( i) = A; = -19.2 MeV, 
= Ai = 39.8 MeV. 



(29) 



with 



|lS2p(| )) - -0.997] J s = |) - 0.080| J s = ±), (30) 
|l52p(|)> - 0.0801 J a = |> - 0.997| J s = \), 
For J = i, the mixing matrix has the form 
/ -23.8 -15.6 ^ 



-15.6 14.2 



with the eigen-vectors given by 



MeV, 



(31) 



|152p(i)> = -0.941|J s =^)-0.338|J s =i),(32) 
|lS2p(i)> = 0.338| J s - |> - 0.941| J, = \), 



and the eigen-values equal 

A' ( ^ } = Af, = -29.5 MeV, 
A^) = A 2 = 19.8 MeV. 



(33) 



We straightforwardly check, that the difference between 
the wave functions as caused by the different masses of 
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TABLE I: The spectrum of diquark levels without the spin-dependent splittings: masses and mean-squared radii. 



Diquark bb 


diquark level 


mass (GeV) 


(r 2 ) 1/2 (fm) 


diquark level 


mass (GeV) 


(r 2 ) 1/2 (fm) 


IS 


9.72 


0.33 


2P 


9.93 


0.54 


2S 


10.01 


0.69 


3P 


10.13 


0.87 


3S 


10.19 


0.99 


4P 


10.30 


1.14 


4S 


10.35 


1.26 


5P 


10.44 


1.39 


5S 


10.49 


1.50 


6P 


10.56 


1.62 


3D 


10.07 


0.72 


4D 


10.24 


1.02 


5D 


10.38 


1.28 


6D 


10.51 


1.52 


4F 


10.18 


0.87 


5F 


10.33 


1.14 


6F 


10.46 


1.40 


5G 


10.27 


1.01 


6G 


10.41 


1.28 


6M 


10.36 


1.15 


Diquark be 


IS 


6.45 


0.48 


3P 


6.91 


1.17 


2S 


6.77 


0.95 


4P 


7.11 


1.52 


3S 


6.99 


1.34 


3D 


6.82 


0.97 


2P 


6.67 


0.75 


4D 


7.03 


1.35 


4F 


6.95 


1.17 


5F 


7.14 


1.53 


5G 


7.07 


1.35 


6H 


7.17 


1.51 


Diquark cc 


IS 


3.13 


0.58 


3P 


3.62 


1.37 


2S 


3.47 


1.12 


4P 


3.85 


1.78 


3S 


3.72 


1.57 


3D 


3.52 


1.14 


2P 


3.35 


0.88 


4D 


3.76 


1.58 



TABLE II: The characteristics of radial wave function for the diquarks: i? d(ns) (0) (GeV 3/2 ), R d{np) (0) (GeV 5/2 ). 



Diquark bb 


Diquark bb 


Diquark be 




Diquark 


cc 


nL 


Rd(ns)(0) 


nL 


R'd(np)(Q) 


nL 


Rd(ns){0) 


nL 


R'd(np)(Q) 


nL 


Rd(ns)(0) 


nL 


R'd(np)(Q) 


nL 


Rd(ns) (0) 


nL 


R'd(np)(fy 


IS 


1.345 


2P 


0.479 


3S 


0.782 


4P 


0.585 


IS 


0.722 


2P 


0.200 


IS 


0.523 


2P 


0.102 


2S 


1.028 


3P 


0.539 


4S 


0.681 


5P 


0.343 


2S 


0.597 


3P 


0.330 


2S 


0.424 


3P 


0.155 



diquark subsystem, is unessential, so that for the 2S2p- 
level the splittings are very close to the values calculated 
above. 

The splitting of diquark, A^ Jd \ is numerically small: 
|A( Jd )| < 10 MeV. Such the corrections are unessential 
up to the current accuracy of method (SM « 30 — 40 
MeV). They can be neglected for the excitations, whose 
sizes are less than the distance to the strange quark, i.e. 
for the diquarks with small values of principal quantum 
number. 

The mass spectra of Qbb and "Ebb baryons are compared 
in Table [II, wherein we restrict ourselves by the presen- 
tation of S-, P- and D-wave levels. 

The most reliable predictions are the masses of baryons 
151s (J p = 3/2+, 1/2+), 2Pls (J p = 3/2", 1/2") 
and 3Pls (J p = 7/2+, . . .1/2+). The 2Pls-level is qua- 
sistable. In the S&6 family the transition into the ground 
state requires the instantaneous change of both the or- 



bital momentum and the summed spin of quarks inside 
the diquark. Therefore, when the splitting between 2Pls 
and lS2p, AE ~ Aqcd, is not small, their mixing is 
suppressed as SV/AE ~ 1 ^~KE ^ Since the 
admixture of \S2p in the 2Pls-state is low, the 2Pls- 
levels are quasistable, i.e. their hadronic transitions into 
the ground state with the emission of 7r-mesons are sup- 
pressed as we have derived, though an additional sup- 
pression is given by a small value of phase space. In 
contract to the Sbb family, in the flbb system the tran- 
sition of 2Pls to the ground level with the emission of 
pion is forbidden because of the conservation of flavor in 
the strong interactions, while the emission of kaon with 
the transition in to the ground state of S^f, is forbidden 
by nonsufficient phase space. An alternative possibility 
is the transition under the emission of two pions in the 
singlet of the isospin group, that is the open channel for 
the 251s levels. 
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TABLE III: The mass spectrum of baryons in comparison with Em, one. 



(ndLciniLi) J p 


M[Q, bb ], GeV 


M[E bb ], GeV [If] 


(ndLdniLi) J p 


M[Q, bb ], GeV 


M[E bb ], GeV [L6| 


(IS ls)l/2+ 


10.210 


10.093 


(3P Is) 1/2" 


10.617 


10.493 


(IS ls)3/2+ 


10.257 


10.113 


(3D ls)5/2'+ 




10.497 


(2P Is) 1/2" 


10.416 


10.310 


(3D ls)7/2+ 


10.627 


10.510 


(2P ls)3/2~ 


10.462 


10.343 


(3P ls)3/2" 


10.663 


10.533 


(2S Is) 1/2+ 


10.493 


10.373 


(IS 2p)l/2" 


10.651 


10.541 


(2S ls)3/2+ 


10.540 


10.413 


(1S 2p)3/2" 


10.661 


10.567 


(3D ls)5/2+ 




10.416 


(1S 2p)l/2'- 


10.700 


10.578 


(3D ls)3/2'+ 




10.430 


(1S 2p)5/2" 


10.670 


10.580 


(3D Is) 1/2+ 


10.617 


10.463 


(1S 2p)3/2'" 


10.720 


10.581 


(3D ls)3/2+ 




10.483 


(3S Is) 1/2+ 


10.682 


10.563 



As for the higher excitations, the 3Pls-states are close 
to the 15'2p-levels with J p = 3/2", 1/2", so that the op- 
erators changing both the orbital momentum of diquark 
and its spin, can lead to the essential mixing with an 
amplitude SV nn > / AE nn > ~ 1, despite of suppression by 
the inverse heavy quark mass and small size of diquark. 
The mixing slightly shifts the masses of states. The most 
important effect is a large admixture of lS2p in 3Pls. It 
makes the state to be unstable because of the transition 
into the ground lSls-state with the emission of gluon 
(the El-transition). This transition leads to decays with 
the emission of kaons 2 . 

The level \S2p J p — 5/2" has the definite quan- 
tum numbers of diquark and light quark motion, because 
there are no levels with the same values of J p in its vicin- 
ity. However, its width of transition into the ground state 
of Sf,b and kaon is not suppressed and seems to be large, 
T ~ 100 MeV. 

One could expect the transitions of 



hi, 



9. 



t,h 



^bb ( - ) K in S - wave, 



-hi, | - ) K in D — wave, 



■=•"» I 2 ) K in D — wave, 



-/,/, | - | if in S — wave. 



The D-wave transitions are suppressed by the ratio of low 
recoil momentum to the mass of baryon. 



2 Remember, that the Sqq/ -baryons are the iso-dublets, while the 



The width of lSls state with J p = 3/2 + is com- 
pletely determined by the radiative electromagnetic Ml- 
transition into the basic J p = l/2 + state. 



B. Q.c C baryons 

The calculation procedure described above leads to the 
results for the doubly charmed baryons as presented be- 
low. 

For lS2p, the splitting is equal to 



A ( i) = 5.1 MeV. 



(34) 



For J = I, the mixing is determined by the matrix 



MeV, 



/ -18.5 -14.8 
I -14.8 42.9 

so that the eigen-vectors 

|lS2p(|')> = -0.975| J s = |) 
|l52p(|)) - 0.223| J s = |) - 0.975| J s 



(35) 



0.223| J s = -),(36) 



have the eigen-values 

A /( ^ = A; = -21.9 MeV, 
= Ai = 46.3 MeV. 

For J = i, the mixing matrix equals 



(37) 



-32.6 -47.6 
-47.6 -31.4 



MeV, 



where the vectors 



-0.712|J S 



0.703|J S 



(38) 



(39) 



ones are the iso-sinelets. 



\IS2 P { 1 -)) = 0.703| J a = |> - 0.712| J s = \), 
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TABLE IV: The mass spectrum of £l cc baryons in comparison with S cc one. 



(udLdniLi) J p 


M[fi cc ], GeV 


M[5 CC ], GeV [If] 


{ndLdfiiLi) J p 


Mpco], GeV 


M[S CC ], GeV [16| 


(IS Is) 1/2+ 


3.594 


3.478 


(3P Is) 1/2" 


4.073 


3.972 


(IS ls)3/2+ 


3.730 


3.61 


(3D ls)3/2'+ 




4.007 


(2P ls)l/2" 


3.812 


3.702 


(IS 2p)3/2' _ 


4.102 


4.034 


(3D ls)5/2 + 




3.781 


(IS 2p)3/2~ 


4.176 


4.039 


(2S Is) 1/2+ 


3.925 


3.812 


(IS 2p)5/2" 


4.134 


4.047 


(3D ls)3/2+ 




3.83 


(3D ls)5/2'+ 




4.05 


(2P ls)3/2~ 


3.949 


3.834 


(IS 2p)l/2'" 


4.145 


4.052 


(3D Is) 1/2+ 


3.973 


3.875 


(3S Is) 1/2+ 


4.172 


4.072 


(IS 2p)l/2" 


4.050 


3.927 


(3D ls)7/2+ 


4.204 


4.089 


(2S ls)3/2+ 


4.064 


3.944 


(3P ls)3/2" 


4.213 


4.104 



have the eigen-values 



A 2 



-79.6 MeV, 
15.6 MeV. 



(40) 



For the 15-, 25- and 3>!5-wave levels of diquark, the 
shifts of vector states are equal to 

A(1S) = 6.4 MeV, 
A(2S) = 4.7 MeV. 
A(35) = 4.2 MeV. 

The mass spectra of the f2 cc and S cc baryons are pre- 



sented in Table IV 



C. flbc baryons 

As we have already mentioned in the Introduction, the 
heavy diquark composed of the quarks of different fla- 
vors, turns out to be unstable under the emission of soft 
gluons. We suppose that the calculations of masses for 
the excited flbc baryons are not so justified in the given 
scheme. Therefore, we present only the result for the 
lowest states with J p — l/2 + 

M n > bc = 6.97 GeV, M Qbc = 6.93 GeV, 

and with J p = 3/2+ 

Mn tc = 7.00 GeV, 

whereas for the vector diquark we have assumed that 
the spin-dependent splitting due to the interaction with 
the strange quark is determined by the standard contact 
coupling of magnetic moments for the point-like systems. 



IV. CONCLUSION 

In this paper we have evaluated the spectroscopic char- 
acteristics of baryons containing two heavy quarks and 



the single strange quark, in the framework of potential 
model. The calculations have been based on the assump- 
tion of string-like structure of doubly heavy baryon, when 
the diquark of small size interacts with the strange quark 
in the limit of quark-diquark factorization in the wave 
functions. We have explored the QCD-motivated model 
of static potential Q , that takes into account two known 
asymptotic regimes at small and large distances. The 
first limit is the asymptotic freedom up to three-loop run- 
ning of coupling constant consistent with the measure- 
ments of a s at large virtualities. The second regime is the 
linearly raising confining potential. The spin-dependent 
corrections have been taken into account. The region 
of factorization applicability as well as the uncertainties 
have been discussed. 

Below the threshold of decay into the heavy baryon and 
heavy strange meson, we have found the system of excited 
bound states, which are quasistable under the hadronic 
transitions into the ground state. In the baryonic sys- 
tems with two heavy quarks and the strange quark, the 
quasistability of diquark excitations is provided by the 
absence of transitions with the emission of both a single 
kaon and a single pion. These transitions are forbidden 
because of small splitting between the levels and the con- 
servation of the iso-spin and strangeness. Further studies 
on the electromagnetic and hadronic transitions between 
the states of doubly heavy baryons are of interest. 

In conclusion we compare the results obtained in the 
present paper with the estimates in potential models and 
in lattice simulations as shown in Table 

The quark-diquark factorization in calculating the 
masses of ground states for the baryon systems with two 
heavy quarks was also considered in ref. ]3^] , where the 
quasi-potential approach 38 was explored. There is a 
numerical difference in the choice of heavy quark masses, 
that leads to that in |37) the mass of doubly charmed 
diquark, say, about 100 MeV greater that the mass used 
in the above calculations. This difference determines the 
discrepancy of estimates for the masses of ground states 
presented in this paper and in |37| . We believe that this 
deviation between the quark masses is caused by the use 
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TABLE V: The masses of ground states M (in GeV) for 
the baryons with two heavy quarks calculated in various ap- 
proaches (* denotes the results of authors in this work). The 
accuracy of predictions under the variation of model param- 
eters is about 30-50 MeV. The systematic uncertainties are 
discussed in the text. 



utxi y u±± 


* 


H7I 




t in 1 
1 ±u ^ 


411 


[121 


4"1 


145 


I4fil 


i — <cc 


3.48 


3.66 


3.74 


3.66 


3.61 


3.65 


3.71 


3.62 


3.57 


1 — 'cc 


3.61 


3.81 


3.86 


3.74 


3.68 


3.73 


3.79 


3.73 


3.63 




3.59 


3.76 


3.76 


3.74 


3.71 


3.75 


3.89 


3.78 


3.69 




3.73 


3.89 


3.90 


3.82 


3.76 


3.83 


3.91 


3.87 


3.75 


1 — 'bb 


10.09 


10.23 


10.30 


10.34 






10.43 


10.20 




— 6S> 


10.11 


10.28 


10.34 


10.37 






10.48 


10.24 




ilbb 


10.21 


10.32 


10.34 


10.37 






10.59 


10.36 




fi» 


10.26 


10.36 


10.38 


10.40 






10.62 


10.39 






6.82 


6.95 


7.01 


7.04 






7.08 


6.93 


6.84 


m/ 

— ct 


6.85 


7.00 


7.07 


6.99 






7.10 


6.96 


6.83 


"ct 


6.90 


7.02 


7.10 


7.06 






7.13 


6.98 


6.88 


Qcb 


6.93 


7.05 


7.05 


7.09 






7.23 


7.09 


6.95 


Q' cb 


6.97 


7.09 


7.11 


7.06 






7.24 


7.12 


6.94 


^cb 


7.00 


7.11 


7.13 


7.12 






7.27 


7.13 


6.98 



of Cornell potential with the constant value of effective 
coulomb exchange coupling in contrast to the above con- 
sideration with the running coupling constant, that can- 
cels the uncertainty in the arbitrary additive shift of en- 
ergy. Furthermore, in the potential approach the masses 
of heavy quarks depend on the mentioned additive shift, 
which adjusted in the phenomemological models by com- 
paring, say, the leptonic constants of heavy quarkonium 
calculated in the model with the values known from ex- 
periments. In the QCD motivated potential such the 
ambiguity of potential because of the additive shift is ab- 
sent, so that the estimates of heavy quark masses have 
less uncertainties. Let us stress that in the Cornell model 
the leptonic constants were calculated by taking into ac- 
count the one-loop corrections caused by the hard gluons. 
This correction is quite essential, in part, for the charmed 
quarks. The two-loop corrections are also important for 
the consideration of leptonic constants in the potential 
approach Moreover, in j37[ the constituent mass of 
light quark is posed with no correlation with the normal- 
ization of potential, while we put the constituent mass 
to be a part of nonperturbative energy in the potential. 
This can lead to an additional deviation between the es- 
timates of baryon masses about 50 MeV. Taking into ac- 
count the above notes on the systematic differences, we 
can claim that the estimates of ground states masses for 
the baryons with two heavy quarks in |57| agree with the 
values obtained in the presented approach (see Table |v|). 

In ref. |^] , following in the framework of quasi- 
potential approach, the analysis of spin-dependent rela- 
tivistic corrections was performed so that the overesti- 
mated, to our opinion, value of heavy diquark mass from 



was used. Unfortunately, there is an evident mis- 
take in the description of calculations in because 
both the parameter giving the relative contribution of 
scalar and vector part in the potential and the anomalous 
chromo-magnetic moment of heavy quark are denoted by 
the same symbol, that leads to numerical errors, since in 
it was shown that these quantities have different val- 
ues. This mistake enlarges the uncertainty about 100 
MeV into the estimates of |59|, so that we can consider 
that the results of |}{| do not contradict with the pre- 
sented description (see Table [v]). 

The estimates based on the hypothesis of pair interac- 
tions were presented in ref. ]40| , so that in the light of 
discussion given in the Introduction the difference about 
200-300 MeV, that follows from values in Table [v|, is 
not amazing. This deviation is, in general, related with 
the different character of interquark forces in the doubly 
heavy baryon, though the uncertainty in the heavy quark 
masses is also important. 

In ref. |yj simple speculations based on the HQET 
with the heavy diquark were explored, so that the esti- 
mates depend on the supposed mass of diquark composed 
of two heavy quarks. In this way, if we neglected the 
binding energy in the diquark, that is evidently related 
with the choice of heavy quark masses, then we got the 
estimates of ground state masses shown in Table [v|. 

Next, in the analysis given in |Q was modified on 
the basis of interpolation formulae for the mass of ground 
state with account for the dependence of spin forces on 
both the wave functions and the effective coupling con- 
stant, which were changed with the quark contents of 
hadrons. In this way, the parameter of energy shift en- 
ters the fitting function, so that this parameter essen- 
tially changes under the transition from the description 
of mesons to baryons: 6m ~ 80 MeV — > 6 b ~ 210 
MeV. This shift of energy provides a good agreement of 
fitting with the mass values for the mesons and baryons 
observed experimentally. However, if we suggest that the 
doubly heavy baryon is similar with the meson contain- 
ing the local heavy source in the picture of strong inter- 
actions, then we should use the energy shift prescribed 
to the heavy mesons but the heavy baryons, wherein the 
presence of system with two light quarks leads to the es- 
sential difference in the calculation of bound state masses, 
hence, to the energy shift different from the mesonic one. 
Such the substitution of parameters would lead to a more 
good agreement between the results of (4^] (see Table ^) 
and the values obtained in this work. 

Recently the analysis presented in |37j was modified 
in in order to take into account the running of QCD 
coupling in the static potential in the framework of quasi- 
potential approach. The authors of |^5| claim that the 
difference between the estimates of masses for doubly 
heavy baryons in Q and Q accumulates three sources. 
The first is the choice of heavy quark masses, as we have 
mentioned above, that gives the shift about 50 MeV by 
the estimates in p5[ . The second source is the relativis- 
tic dispersion law of free light or strange quarks. The 
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third is the spin-independent relativistic corrections for 
the orbitally excited states. Nevertheless, we emphasize 
that the estimates obtained in these two approaches are 
in a good agreement within the limits of systematic un- 
certainties about 70 MeV. 

Finally, the lattice simulations based on the lagrangian 
of NRQCD are presented in Table 0, too Q. One can 
see that the lattice results show approximately twice re- 
duction of spin-dependent splitting of the ground level, 
which position agrees with the estimates obtained in the 
potential models. 

Summarizing, we can claim that, first of all, in the 
framework of potential approach in the calculations of 
masses for the doubly heavy baryons the dominant un- 
certainty is caused by the choice of heavy quark masses, 
so that due to the adjustment on the systems with heavy 
quarks, the analysis presented in the QCD-motivated 
model of potential with the running coupling constant at 
short distances and the linear nonperturbative term con- 



fining quarks at large distances, gives the most reliable 
predictions. Further, the calculations in the framework 
of NRQCD sum rules jL7) gave the results for the ground 
state masses with no account of spin-dependent forces, 
so that the sum rule estimates are in a good agreement 
with the values obtained in the potential models. 
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